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Abstract 

A quantum field theoretic treatment of inclusive deep-inelastic diffractive scattering is 
given. The process can be described in the general framework of non-forward scatter- 
ing processes using the light-cone expansion in the generalized Bjorken region. Evolution 
equations of the diffractive hadronic matrix elements are derived at the level of the twist- 
2 contributions and are compared to those of inclusive deep-inelastic forward scattering 
(DIS). The diffractive parton densities are obtained as projections of two-variable par- 
ton distributions. We also comment on the higher twist contributions in the light-cone 
expansion. 



1 Introduction 



Deep inelastic diffractive lepton-nucleon scattering was observed at the electron-proton collider 
HERA some years ago [I]]. This process is measured in detail by now || and the structure function 
F®(x, Q 2 ) was extracted^ The diffractive events are characterized by a rapidity gap between the 
final state nucleon and the set of the diffractively produced hadrons due to a color-neutral 
exchange. The (semi-inclusive) structure functions emerging in the diffractive scattering cross 
section show the same scaling violations as the structure functions in deep inelastic scattering 
within the current experimental resolution H|. This is a remarkable fact, which should also be 
understood with the help of field-theoretic methods within Quantum Chromodynamics (QCD). 
In the description of diffractive ep scattering one has to clearly distinguish || the case of hard 
scattering, i.e. large Q 2 , from that of softer hadronic interactions, cf. also ||. A perturbative 
description of the scaling violations can only by hoped for in the former case, to which we limit 
our considerations in this paper. 

There is a vast amount of varying descriptions of the underlying dynamics of the diffractive 
scattering process^, ranging from intuitive phenomenological models || to non-perturbative semi- 
classical descriptions being able to describe and to parameterize O] the existing data 
differentially. 

In this paper we describe the process of inclusive deep-inelastic diffractive scattering, being 
a non-forward process, at large space-like momentum transfer using the general light-cone 
expansion for the non-forward case, cf. [12-14]. As will be shown, the scaling violations of the 
diffractive structure functions can be described perturbatively in this region. We show that the 
diffractive parton densities can be derived as a special projection of two-variable distribution 
functions f A (z + , where z + and z_ are light-cone momentum fractions; for other projections 
in similar cases related to different observables see e.g. [13-15]. In particular we have no need 
to refer to the specific mechanism of the non-perturbative color-singlet exchange between the 
proton and the set of diffractively produced hadrons and consider instead the expectation values 
of operators of a given twist and their anomalous dimensions. In this way we can generalize the 
analysis to operators of higher twist. This formulation is very appropriate for potential later 
studies of the corresponding operator matrix elements with the help of lattice techniques similar 
to the case of deep-inelastic scattering [16| calculating their Mellin-moments. 

The paper is organized as follows. In section 2 we consider the Lorentz structure of the 
differential cross section for inclusive deep-inelastic diffractive scattering. The short-distance 
structure of the matrix element is derived in section 3 using the (non-local) light-cone expansion 
in the generalized Bjorken region. In section 4 we derive the anomalous dimensions for the case 
of twist-2, discuss the implications for higher twist operators and compare the case of diffractive 
scattering to that of deep-inelastic scattering, and section 5 contains the conclusions. 



2 Lorentz Structure 

The process of deep-inelastic diffractive scattering is described by the diagram Figure 1. The 

lr The measurement of the longitudinal diffractive structure function F£(x,Q 2 ) has not yet been possible. For 
the DIS structure function cf. ||. 

2 For surveys and a recent account see |7j [|. 
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differential scattering cross section for single-photon exchange is given by 



2(s - M 2 ) 4 



spins 



(1) 



Here s = (pi + I) 2 is the cms energy of the process squared and M denotes the nucleon mass. 

V 
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Figure 1: The virtual photon-hadron amplitude for diffractive ep scattering 

The phase space dPS^ depends on five variables^ since one final state mass varies. They can 
be chosen as Bjorken x = Q 2 /(W 2 + Q 2 — M 2 ), the photon virtuality Q 2 = —q 2 , t = (pi — p 2 ) 2 , 
a variable describing the non-forwardness w.r.t. the incoming proton direction, 



x ¥ 



2r) 



Q 2 + M 2 x -t 

I-77 Q 2 + W 2 - M 2 



> x , 



(2) 



demanding M\ > t and where 



<1-(P2-Pi) 
V = — e 



2 - x 



Q-(P2+Pl) 

and $ the angle between the lepton plane p x x I and the hadron plane p l x p 2 , 

(Pi x 0-(Pi x Pa) 



(3) 



cos $ 



|Pl X '||Pl x P2I 



(4) 



^ 2 — G°i + l) 2 an d M£ = (pi + q — P2) 2 denote the hadronic mass squared and the square of 
the diffractive mass, respectively. 

Since the leptonic tensor L^ u is symmetric for unpolarized scattering and the electromag- 
netic current is conserved the hadronic tensor, which is generally composed out of the tensors 
9pv, {q.hiPiij.iP2h) ® (<lv,Piv,P2v), consist out of four structure functions onlyQ 
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Q 
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Piv - q v ^r 
T 
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(5) 



3 For the commonly used notation of four variables see e.g. |Tj 
4 For a construction method, see e.g. J18[. 
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with 



Wi = Wi(x,Q 2 ,x F ,t) . 



(6) 



The size of the rapidity-gap [19[ || is of the order Ay ~ ln(l/xp) and large in the kinematic 
domain of HERA. 

If target masses can be neglected, M 2 ~ 0, and t is considered to be very small compared to 
all other invariants the in- and outgoing proton four-momenta become proportional, p 2 = zp\. 
The variable z is then related to Xf and r] by 



z = 1- x ¥ = - . 

1—7] 

In this limit the hadronic tensor is described by two structure functions 



W, 



fJ,U 



Pip - J [Piu - g v --^ 



Pi-q\ W2_ 

M 2 ' 



with 



(7) 



W 2 = W 3 + (1 - x P )W 5 + (1 



'Wa . 



(9) 



Due to the dependence on xp or rj, Eq. (0), the process is non-forward although the algebraic 
structure of the hadronic tensor is the same as in the forward case. Finally the generalized 
Bjorken-limit is carried out, 



2p 1 .q = 2Mv — > oo, p 2 .q — > oo, Q 2 — > oo with x and 



fixed 



(10) 



which leads to MWi(x,Q 2 ,x v ) — > F 1 (x,Q 2 ,x P ) and uW 2 (x, Q 2 , x P ) — > F 2 (x, Q 2 ,x P ). Note that 
Q 2 3> t,M 2 always holds in this region. In the limit p 2 = zp\ and M 2 ,t = the $— integral 
becomes trivial. The scattering cross section reads 



O'diSv 



2ira z 



dxdQ 2 dxf> xQ 4 
with y = q.pi/l.pi. 



y 2 2xF? {S \x, Q 2 , x r ) + 2(1 - y)F? w (x, Q 



D(3), 
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3 The Compton Amplitude 



The renormalized and time-ordered product of two electromagnetic currents is given by 



iRT 



—e 



J - (§) J » 



X 



2 

RT 



2tt 2 (x 2 — ie) 2 
x denotes a light-like vector corresponding to x 

c 



s 



5(12) 



x = x + 



e 



a/ x.( 2 — x 2 ( 2 — x.( 



(13) 
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and ( is a subsidiary vector. Following Refs. [14, 20] the operator T^ u can be expressed in terms 
of a vector and an axial-vector operator by 



x _x\ I X _x 

2' 2} £ ^ u ° ">\2' 2 



t7r*{x* — it) 
where 

Safj,Xu 9afj,9\u 9Xfi9au 9 fiv9Xa • 

The bilocal light-ray operators are 



O c 



x x 
2'~2 
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(14) 

(15) 

(16) 
(17) 



The general operator is now to be related to the diffractive scattering cross section. This is 
possible upon applying Mueller's generalized optical theorem |2]J (Figure 2), which moves the 
final state proton into an initial state anti-proton. 
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Figure 2: A. Mueller's optical theorem. 

Consequently, the Compton amplitude is obtained as the expectation value 

T^(x) = (p h -p^T^px, -p 2 ) , (18) 

which is forward w.r.t. to the direction defined by (pi, — p 2 \. The quantity (pi, — p 2 \, in the 
strict sense, is not a single physical state of the Hilbert-space of states, since this would imply 
t = (P2 ~ Pi) 2 > 0. The analytic continuation is such, that t is kept space-like, which is 
necessary for the non-perturbative behavior of the matrix element. The twist-2 contributions 
to the expectation values of the operators (|l(||17]) are obtained 



($i,-P2\Of^{K+x,K-x)\px,-p 2 ) = / dXd^ipt, -p 2 \Of B) (\n + x, Xn-x^pt, -p 2 ) 

J 



(19) 



as partial derivatives of the expectation values of 

O a (k + x, k_x) = x a O A (n + x, k_x) (20) 
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the corresponding scalar and pseudo-scalar operators. The index A = q,G labels the quark- or 
gluon operators, cf. f[4" |. In the unpolarized case only the vector operators contribute. The scalar 
twist-2 quark operator matrix element has the representation^ due to the overall symmetry in x 

(p x , -p^O^K+x, K-x)\p u -p 2 ) = xp^Dz e~ iK - xp * /«(*+, zJ) + xn^Dz e~ iK ~ xp * tf(z + , (21) 

and is independent of k + . Here we assume that all the trace-terms have been subtracted, see 



rjjj. The action of (pi, — p 2 1 onto the scalar operator yield projections onto two directions p_ 
and 7r_ = p + — p~/r). f A {z + ,zJ) and f£(z+,z-) denote the corresponding scalar two-variable 
distribution amplitudes and the measure Dz is 

Dz = dz + dz-9(l + z + + z-)9(l + z + - Z-)9(l - z + + zJ)6(l - z+ - zJ) . 

Here, we decomposed the vector p z as 

p z = P-Z- + p+Z + = + TI-Z+ , 

with zi t 2 momentum fractions along px,2 and p± = p 2 ±Pi, z± = (z 2 ± zi)/2 and 



(22) 



(23) 



$ = z_ H — z + . 
V 



(24) 



Note that, q.n_ = 0. In the approximation M 2 , t ~ 0, in which we work from now on, the vector 
7r_ vanishes. In this limit only the first term contributes to the matrix element Eq. fl2"T|). 
The Fourier-transform of the Compton amplitude is given by p0| 



, l d^xe^T^x) 



2S p ^ au I Dz 



F(z + ,z.) , (25) 



with Q z = q — Pz/2. The function F(z + , zJ) is related to the distribution function f(z+, zJ) by 

* x d\ „ ( z± Z-' 



The denominators take the form 



1 



Ql + is qp- ($-2/3 + is) 



where 



P 



x 



r 



x ¥ 2q.p_ 

The conservation of the electromagnetic current is easily seen 

9% 1 /(Pi,P2,g) = q u T fll/ (p 1 ,p 2 ,q) = -2p_ ( „ M J DzF( 



(26) 



(27) 



(2* 



z + ,z_) = 



(29) 



3 For parameterizing hadronic matrix elements see e.g. j22| 
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since the symmetry relation, |20 



DzF(z + ,z_) = 
holds. Subsequently we will use the distribution 

F(d, 77) = / DzF(z+, z_)5{<d -z-- z+/7}) ■ 



(30) 



sign(i?)/7) ^ z ^ 

— f(z,v) ■ 

1? z 



(31) 



The distribution function f(z,T]) is related to f(z + ,Z-) by 

,7,(1-*) 



/(*, 77) = / dp 6(1 - p)8(p + l)f(p, z - p/ v ) 

>ri(l+z) 



(32) 



with p = z + /rj. 

The Compton amplitude may be re-written as 



r-Vv ( 1 
/ d$< [-q.p-g^ - 20p_ iM p_,„ + (p_, M g„ + P-,^)] n2 - . 

■/+l/r? [ V« + !£ 



(Ql + 

sign(i?)/rj ^ z ^ 

—f(z,v) ■ 



(33) 



The -^-integrals in Eq. (|33"D can be simplified using the identities 



-!/»? flk / , sign(i?)/r7 ^ ^ 

, 7 „ "v-w+urh —f(. z 't) 



-1/u fctfk-l rri&W/r, dz 

av- r / — tiz.r)) 



- 1/V M .^7>,r|) 



+1/*) 



(1? - 2/3 + ie) 



(34) 



One may further re- write Eq. (|33|) referring to p± instead of p_ 



T tl v{pi,P2,q) = 2 



-l/jj 



~y/xi> H r !/Ip 



(0 - 2/3 + is) 



(35) 



Before we consider the absorptive part of the Compton amplitude we exploit the symmetry 
relation for the unpolarized distribution functions F (zi,z 2 ), Ref. Ml, 



It translates into 



and, cf. Eq. (0), 



F A (z 1 ,z 2 ) = -F A (-z 1 ,-z 2 ) 



F a ($,t 1 ) = -F A (-d,7]) 



(36) 



(37) 



f A {$,v) = -f A (-$,v) ■ 



(38) 



Seeking the form of Eq. (§) , Eq. (|35|) transforms to 



T fiu (p 1 ,p 2 ,q) 



-l/r, 



d$ 



l/r, 



+ 



2x 



pi.q 



q.pi 

f(-#,v) 

d-2(3 + ie d - 2(3 + is 

■l/r. 



Piti - \V\v- q» — 



Pi-Q 



q- 



+ 2 



Q-Pi J l/r, 

Taking the absorptive part one obtains 
1 



d-d ($x P - 2x) 



ft- 2(3 + is 



(39) 



lix 

V Q J Q-Pi V 



(40) 



Pl-Q 



piu - 0"~r) F 2 ( y p,ri,Q 2 



where 



(41) 



F 1 ((3, rj, Q 2 ) = J2e 2 q [f?(P, Q\ x,) + />, Q 2 , a*)] = if (3) (x, 

9=1 

with Nf the number of flavors, choosing the factorization scale /i 2 = Q 2 , and 

F 2 (/?,r/,Q 2 ) = 2xF 1 ((3,r ] ,Q 2 ) = F^ 3 \x,Q 2 ,x P ). (42) 

If we compare the distribution functions fjy with the parton distributions in the deep inelastic 
case, the role of x in the latter case is taken by the variable (3. Both variables have the range 
[0, 1], through which Mellin transforms w.r.t. these variables can be evaluated. On the contrary, 
the support in x of the diffractive structure functions is limited to [0,xp]. Although (3 in the 
diffractive case compares to x in the deep inelastic case, the Callan-Gross relation, at lowest 
order in a s , Eq. (|42|), is given by the factor 2x between the two structure functions. 
The diffractive quark and anti-quark densities are given by 



9=1 



9=1 



-7(2/3, rj, Q 2 ) 
f(-2(3, V ,Q 2 ) 



(43) 



We finally can express the diffractive parton densities in terms of the distribution function 
f(z+, z_) directly 



f u ((3,Q 2 } x P ) 



f (P,Q 2 ,x ¥ ) 



2 — a^p 
xm> — 2x 



dpf{p, 2(3 + p(2-x P )/x P ,Q 2 ) 



xm> J r-2x 



dpf{p, -2(3 + p{2~x F )/x P ,Q 2 



(44) 
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4 Evolution Equations 

The evolution equation of the scalar twist-2 quark and gluon operator Eq. (|19|) read, cf. |L4| .[] 
^ 2 ——0 A ( y K + x,K_x;jj 2 )= [ Dk''J AB (k + , k_, k' + , k'_; n 2 )0 b{k' + x , k'_x; /j 2 ) , (45) 



djj, 

with k± = («2 — k{)/2 and the measure 

Dk = dK + dK-9(l + k + + «_)0(1 + k + - «_)0(1 - k + + «_)0(1 - «+ - «_) . (46) 

We consider the non-forward evolution equations to trace an eventual ^-dependence. Here 
j AB (tz + , K-, k' + , k'_, /i 2 ) is the general non-forward twist-2 singlet evolution matrix. [] The same 
evolution equation applies to the matrix elements (p±, — p 2 \O a (k + , K-)\pi, — p 2 ), which may be 
related to the distribution functions f A ($,rj) byQ 

f A V,l) = f e iK -^{p 1 ,-p 2 \0 A (K + ,K-)\p 1 ,-p 2 )(xp-) 1 - d A , (47) 



2tt 

where d q — 1 and do = 2. Eq. fl47D and the inverse Fourier transform 

(pi,-P2\0 A (K + ,K^)\p ll -p 2 }(xp_) 1 - d *= [ XN dd> tr*-**-*^^) (48) 

Jl/r, 

are used to obtain evolution equations for f A (-d,r}]/i 2 ). To perform the ^.-integral the depen- 
dence on this variable has to be made explicit in the anomalous dimensions ^ AB to all orders. The 
scalar operator matrix elements, Eq. fl2lD , do not depend on the light cone mark k+ due to trans- 
lation invariance on the light-cone, which therefore can be set to zero in Eq. (EDI). Furthermore, 
the anomalous dimension j AB obeys the re-scaling relation, cf. [14| , 



7 AB (k + , k_, k' + , k'_; h 2 ) = a dAB ~i AB (oK +1 <tk_, ctk' + , gk'_) , (49) 
with d^B — 2 + d^ — ds- The anomalous dimension transforms to 

n / d B -d A AB fn i K + K - 2 I f n d B -d A f>AB/ 2\ 

Dk k_ a 7 10,1,—, — ; /i J = / Da n_ A K (ai,a 2 ;/i) 



Def 



d B -d A 



/ cfit(l — w) / <i£ 
Jo io 



x/^(£(i- u ),(i-0(i-«);^) 

= K d _ B ~ dA / dwK (u; /i 2 ) , (50) 
./o 

with a x = [1 - « + k'_)/k_]/2 = £(1 - u), a 2 = [1 - « - k'_)/kJ\/2 = (1 - 0(1 - «)■ The 
/t_-integral yields 

r dK ,v ~ f for A = J B = g,G 

/ e i K _ ip _ { ^ u #>)^ p _ K _y B - dA = qAB( u # = ) d u 5{ud' -■&) for A = q, B = G 

J 271 { 9{u$' -■&)/$ forA = G,B = q 

6 The non-singlet evolution equations are structurally the same as that with the anomalous dimension j qq . 
7 The leading order and next-to-leading order non-singlet and singlet anomalous dimensions were calculated 
in Refs. @ gj| and @. 

8 Note that we identified here r = x.p-/x.p+ with r/ = q.p_/q.p+ which is possible after the Bjorken limit is 
taken. 
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The following evolution equations for the distribution functions f A ($, 77) are obtained 



A /•! P- sign(i?)/7) _ ^AB 

'• ,. 1 , ■■> / ^ / ,,/.-; 1 ;> , . , ; . 

'0 J-d 



V 2 J^fM V, H 2 ) = I du I dtf AB (u$' — ■&) K («; /. 2 )/b(^', r?; /i 2 ) . (51) 



The functions 



# J du AB (u$' - $)K AB (u; fi 2 ) = P AB f^j, A (52) 
are the forward splitting functions, which are independent of 77 resp. xp, leading to 

j p- sign(i?)/r? j q/ / o \ 

^/W)=jf £p»(±,s)M*,*f). (53) 

They act on the momentum fraction d of the diffractive parton densities. Note that the range 
of i?' is yet different from that of the corresponding quantity for deep inelastic scattering. After 
identifying ■& — 2(3 taking the absorptive part, however, one arrives at the twist-2 evolution 
equation 

S±fi(P,x r ;f)=J 1 yjrP2 (^;/) fE(P',x P -^ 2 ) . (54) 

In the case of deep inelastic scattering the corresponding identification for the momentum fraction 
is z = x. The dependence of the diffractive parton densities w.r.t. 77, or xp, is entirely parametric 
and not changed under the evolution, which affects (3. In the case of the twist-2 contributions 
factorization proofs were given |25|P|, leading to the same evolution equations, (53), which are 
confirmed by the present derivation. For phenomenological applications to next-to-leading order, 
see e.g. |f28| . 



We expressed the Compton amplitude with the help of the light-cone expansion at short 
distances and applied this representation to the process of deep-inelastic diffractive scattering 
using Mueller's generalized optical theorem. This representation is not limited to leading twist 
operators but can be extended to all higher twist operators synonymously. The corresponding 
evolution equations for the higher twist hadronic matrix elements, which depend on more mo- 
mentum fractions $j than one, transform analogously to the case of twist-2 being outlined above. 
A central parameter k + = - X^Li Ki ma y be set to zero, and analogous re-scaling relations apply. 
By virtue of this also here forward evolution equations are derived. However, the connection of 
the momentum fractions to the outer kinematic parameters is less trivial than in the case of 
twist-2 due to the structure of the corresponding Wilson coefficients. This is the case also for 
deep inelastic scattering, cf. e.g. 



5 Conclusions 

The differential cross section of unpolarized deep-inelastic ep-diffr active scattering is described 
by four structure functions for pure photon exchange, which depend on the four kinematic vari- 
ables, x, Q 2 ,xp and t. In the limit of vanishing target masses and t — > only two structure 
functions contribute. In the case of hard diffractive scattering the scaling violations of these 

9 They apply correspondingly to the kinematic of fracture functions, see |p6f and for an analysis in <£> 3 theory 
in d = 6 dimensions [^7) . 
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structure functions can be described perturbatively. This is possible in transforming the non- 
forward amplitude of the process via Mueller's optical theorem, by rotating the final-state proton 
into an initial-state anti-proton. The Compton-amplitude is calculated for the hadronic two- 
particle state (pi, — P2I and w.r.t. this state in the forward direction. The diffractive parton 
densities are associated to the general two-variable distribution functions, which describe the 
hadronic matrix element. The scaling variable of the diffractive parton densities, which directly 
compares to the Bjorken variable x in the deep-inelastic case, is /3 — x/x-p. However, the Callan- 
Gross relation takes the usual form F®(x, Q 2 , xp) = 2xFf(x, Q 2 , xp). Due to the transformation 
of the problem, which is made possible applying Mueller's optical theorem, the anomalous di- 
mensions are the same as for forward scattering. We demonstrated this by an explicit calculation 
in the case of the twist-2 operators, however, the same mechanism applies also for higher twist 
operators using the light cone expansion. In the case of the twist-2 contributions the effective 
momentum fraction, $ = z~ + z + /i] may be identified with the variable 2/3 for diffractive scat- 
tering. The dependence of the parton densities on xp is not affected by QCD-evolution, which 
acts on the variable f3. 

Acknowledgment. For discussions we would like to thank J. Bartels, W. Buchmiiller, J. Dain- 
ton, B. Geyer, P. Kroll, and G. Wolf. 
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